Limitations on the Ginzburg criterion for dirty superconductors 
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INTRODUCTION 

The superconducting phase transition has been ver 
well described by the Ginzburg-Landau (GL) theory, [1| 
which is equivalent to the BCS theory [2] for tempera- 
tures T close to the transition temperature T^. However, 
both theories use mean fields, and therefore miss the ef- 
fects of critical fluctuations as T approaches T^. Q A 
phenomenological way to estimate the region of validity 
of these mean-field theories is to calculate a certain ef- 
fect of the fluctuations, and to require that this effect 
be small compared to the corresponding mean-field pre- 
diction. The relative temperature range for which the 
mean-field theory breaks down, Iq = {Tq — T^)/T^, is 
called the "Ginzburg region" . [3| 

Here we discuss dirty superconductors, which contain 
non-magnetic impurities. It is well known that such 
impurities reduce the superconducting coherence length, 
while leaving the transition temperature, as well as the 
density of states, largely unchanged. [1, 0] As a re- 
sult, although the Ginzburg region is still rather small, 
it is significantly larger than for the pure superconduc- 
tors. For T > T^, the superconducting order parameter 
vanishes, and the GL contribution to the specific heat 
vanishes. This specific heat has a discontinuity AC at 
T^. In contrast, the fluctuations in the order param- 
eter generate non-zero contributions to the speciflc 
heat even above T^. One common way to estimate the 
Ginzburg region is then to require that outside of this 
regime one has < AC. [5|] As T approaches T^, 
Cg diverges in c? < 4 dimensions as where 
t = ln(T/Tj « (T - TJ/T^. Keeping only this leading 
divergent term, one finds - (AC/A'')2/(''-4), where 
A is the momentum cutoff, which is the inverse of the 
relevant size of the fiuctuations in space (to be discussed 
below) . 

The above result ignores corrections to the leading di- 
vergent term in Cg . As we discuss below, such corrections 
arise both from quantum fluctuations and from correc- 



tions to the leading 'static' contribution. The aim of the 
present paper is to discuss the effects of these corrections, 
which become crucial as AC becomes small. In order to 
derive these corrections, it is important to obtain the full 
expression for the leading wave- vector and frequency de- 
pendent order-parameter correlation functions, and not 
just the 'static' Ornstein-Zernike (OZ) Q approximation 
x{q) ~ 1/(9^ + ?~^) which is used in the 'standard' GL 
theory (^ is the coherence length). The derivation of this 
full expression is reviewed in Sec. II. Section HI then 
presents the resulting contributions of the fluctuation to 
the specific heat, Cg, including all the corrections, and 
Sec. IV discusses the consequences for the Ginzburg re- 
gion. The results are summarized and discussed in Sec. 
V. 



THE PARTITION FUNCTION 

We begin by reviewing the microscopic derivation of 
the free energy which determines the superconducting 
fluctuations. The Hamiltonian is 



H = y" drn{r) , 



(1) 



with 



nir) 



V4(r)Ho(r)^,(r) 



-l/(r)i/^j(r)i/^|(r)i^^(r)i^^(r) , (2) 

where "0 creates an electron with spin a at r. The 
interaction V{r) depends on the spatial coordinate r, 

Vir) = A(r)/AA(r) , (3) 

where A(r) is the effective local (dimensionless) electronic 
coupling, while A/'(r) is the local density of states per unit 
volume and unit energy. The single-particle part of the 
Hamiltonian ^ reads 



Tin 



-V^/{2m) + u{r) - fi , 



(4) 
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where /i is the chemical potential and the disorder po- 
tential ufr) is modeled by point-like non-magnetic scat- 
terers. Q 

The following calculation has been given, for quasi- 
one dimensional rings, in Refs. |^ and 1^. However, we 
repeat it here for the general d-dimensional case, in or- 
der to highlight the origin of the two Matsubara frequen- 
cies which generate the slow momentum and frequency 
variation of the quadratic coefficient in our effective GL 
theory. The quantum partition function Z is [lol | 



Z= \ P{7A(r,T),^(r,r)}exp[-5] , 



(5) 



where the action S is 

f-/3 



5 = y^r^ dr(^Va(r,r)|-V',(r,T)+H(r,- 



(6) 



and/3 = 1/T (we use K — — 1). Here, the annihilation 
and creation field operators in the Hamiltonian ^ (?/; 
and ■0^) are replaced by the Grassmann variables V'(r, t) 
and ■!/;(r,T), respectively. 

Applying the Hubbard- Stratonovich transformation to 
Eq. ([5]) , and integrating the fermionic part of the action, 
the partition function is cast into the form 



2 = 1 P{A(r,r),A*(r,r)}e- 



(7) 



with the action 

Jo V{v) 



S 



Tr 



{ln[/3g-')}, (8) 



where Q ^ ia the inverse (2x2 matrix) Green function 
at equal positions and imaginary times. 



with 



G- 



A 



A* G: 



being the particle inverse Green function, and 



gj: = -dr 



(9) 



(10) 



(11) 



being the inverse Green function of the holes. 

The integration over the bosonic fields in Eq. (O is 
carried out using a stationary-phase analysis [lC| of the 
action S. At temperatures above the transition temper- 
ature, this amounts to expanding the second term on the 
right-hand side of Eq. ([8|) to second order in A (the 
first-order contribution to the expansion being zero) 



Tr{ln(/3^-i)} 



Tr 



{ln/3 







G7 



drdr' drdr' , /na// /\a*/ n 

— n(r,r',T-T')A(r',T')A*(r,T) , 

(12) 



P2 



where J7 denotes the volume of the system. The first 
term on the right-hand side of Eq. p2|) gives the parti- 
tion function of noninteracting electrons; the second one 
represents the contribution of the superconducting fluc- 
tuations to that function. Its calculation requires the 
correlation 

H(r,r', t-t') = 

-(Gp(r,r',r-r')G,(r',r,r'-r)) , (13) 

where (...) indicates averaging over the impurity config- 
urations (see Ref. for details). For the general case, 
when the material contains several different regions (as 
e.g. in a double layer 9]), this average depends sepa- 
rately on r and on r', and the calculation becomes dif- 
ficult. However, for the purposes of the present discus- 
sion, it suffices to consider the simplest case, in which the 
material is homogeneous, and therefore the attractive in- 
teraction F(r) = V as well as the density of states TV 
are constant in space. Averaging over the impurities re- 
stores homogeneity even in the dirty case, and the spatial 
dependence of H becomes a function of r — r'. Hence, 

Tr{ln(/3^-i)}[_^^ =^5]H(q,^)|A(q,^)P , (14) 

V q 

where 

n(q,i') = ^ ^{G{^^ +q,P2 +q,a; + jy) 



P1.P2 ^ 



xG(-pi,-p2,-c^)) , (15) 



and both Green functions are the particle ones, [8| i.e., 
G = Gp. We use the notations w = w„ = 7rT(2n -I- 1) 
for the fermionic Matsubara frequencies, and v = = 
lnTm for the bosonic ones (n and m are integers) . Since 
the phonon-mediated electron-electron attractive inter- 
action is limited to energies within the Debye frequency 
ujy) from the Fermi energy, both |w| and |a;-|-j^| are bound 
by LO^. 

Inserting these results into the expression for the action 
[see Eq. ([H])] , the Gaussian fluctuation-induced partition 
function, 2fl^2, takes the form 



Zfi 2 = y 2?{A(q, ^.),A*(q,^.)}^ 



-52 



with 



52 = 5] 5] A* (q, z.) - n(q, v)) A(q, v) 



(16) 



(17) 



The function H(q, z^), Eq. (|T5|) . is calculated by ex- 
tending the method employed in Refs. |llj and jl2 to 
include the dependence on v and on the wave-vector 
q. The calculation is valid in the diffusive limit, in 
which (D II'kT')^!'^ is much larger than the mean- free 
path I = VpT of the relevant metal, where D = P /{dr) 



3 



is the diffusion coefficient, Vp is the Fermi velocity and r 
is the mean free time associated with the scattering from 
the non-magnetic impurities. This is equivalent to the 
requirement that 2ttTt <Si 1. Alternatively, we require 
that T be the shortest time in the problem, or that all 
the energies (including uijj and T) be smaller than h/r. 



Following the derivation given in Ref. 
the response function 11 in the form 



11 



n(q,:/) = ^X(q,j., 



we present 



(18) 



The function i^(q, z/, w) obeys a diffusion equation (with 
the diffusion constant D) due to the scattering by the 
non-magnetic impurities. It also obeys a Dyson equation, 
which yields 



2TrNQ[Lu{L0 + ly)] 
\2uj + + Dq^ 



(19) 



where TV = ^lAf{0) is the density of states per unit en- 
ergy at the Fermi energy. The appearance of two fre- 
quencies in the denominator of Eq. (|19|) is crucial for 
the discussion below. These two frequencies result from 
the appearance of the two Green functions in Eq. ([T5|) . 
Finally, we find 



where 

a(q,z.,r) = i-^n(q,z.,r) . (21) 

Since the Debye frequency oj^ serves as a cutoff on the 
fermionic Matsubara frequencies uj = 7rT(2n +1) in Eq. 
dH]), one finds 



^n(q,.,T) 



Dq' 



L»q2 + 2w£ 



AttT 



(22) 



where ^ is the digamma function. 

Equation (I^IH) represents the first term in an effective 
GL-like expansion of the free energy density in powers 
of the order parameters A(q, ;/). The original GL the- 
ory [H was phenomenological, and the coefficient a was 
written as a — a'{T — T^) + cq^, which is equivalent 
to the Ornstein-Zernike approximation. This ignored 
the quantum fluctuations, and was presumed valid very 
close to and for long wave-lengths. An extension 
which includes quantum fluctuations had a(q, j/, T) — 



a'{T - TJ -f c|q|2 + e\v\"' /\q\'^ . [l3l Such forms (with 
m — \ and m' — 0) can also be obtained from the micro- 
scopic expression (|2ip , if one expands that expression to 
the lowest order in {T — T^), in q^ and in v. However, it 



is clear that an expansion oi^ mv = 27rTm is not justi- 
fied for any m 7^ 0, and that an expansion in powers of q 
may also be allowed only in a limited range of momenta. 
Much of the literature restricts itself to the 'static' limit, 
i.e. v — Q, and to the expansion of a up to order q^. 
5] Our purpose here is to investigate the full expression 
pTjl . and find deviations from the simpler GL theory. 

Within such an effective Ginzburg-Landau theory, the 
phase transition occurs when the coefficient a(q, T) 
first vanishes as the temperature T is lowered. Using 
Eq. (|22p , this happens for the fluctuation- free "classical" , 
Landau or mean-field limit, q = = 0. For ^ T, we 
use the asymptotic limit ^{z) ^ Inz for large z and find 
the transition temperature 

27b, 



Tr 



-1/A 



(23) 



where 7^ — e^*t^/^l/4 is the Euler constant. This result 
is consistent with the usual BCS one. 

For ijjj^ ^ T (i.e. for small A), we replace 1/A in Eq. 
([2T|) by Eq. For T close to T^, we also denote 

t = \n{T /T^), and replace T by T^, in the denominator of 
Dc^ /{AttT). Substituting also v = 2nTm we find 



a(q, i',T) = t + * 



■ 1 + Iml 



D 





T- 








.2. 



(24) 



This is the expression we shall use in the following calcu- 
lations. In some expressions below we change notation, 
a(q,i/,r) ^ a(q,m,r). 

In the "pure" Landau theory, one keeps only the "clas- 
sical" term, with Aq = A(q = 0,1/ = 0), and one adds 
the quartic term in the free energy density, of order | Aq j"*, 
so that the Landau free energy density has the form 



TlnZ" =r5«A/'(ao|AoP 



1 



(25) 



where Z'^ contains only Aq and (for large ujjj) ag = 
a(0,0,r) = A-i - n(0,0,r)/(/3iV) « t. Microscopic 
calculations yield b = 7C{S)/{8tt'^T'^) = b^/T'^, with 
6q = 0.1. [5I, For metals, we write the electron energy 
as -E = pv/2, and therefore the density of states (per unit 
volume and unit energy) at the Fermi level is given by 
N{Q) = S^f^/vp = Syp/{2Ep), where Ep and Pp 
are the Fermi energy and momentum, 5^ ~ A^/(27r)'*, 
and = 27r''/Vr(d/2) is the area of the unit sphere in 
d dimensions (F is the gamma function) . Minimizing 
with respect to Aq, one finds a non-zero Aq below and 
a jump in the specific heat (per unit volume), from zero 
to 

AC^"^) = A/'(0)/[6Tj - 10A/-(0)r, = bS^pTjEp . 

(26) 

However, the fluctuations at non-zero wave-vectors and 
frequencies give important contributions to the specific 
heat even above T^,, as we discuss next. 



4 



SPECIFIC HEAT DUE TO FLUCTUATIONS 



Substituting Eq. ([20)) into Eq. ([16)) yields a Gaussian 
integral, with the result 



fl,2 



nn 

q V 



1 



.(q,i.,T) 



(27) 



where unimportant multiplicative factors have been 
omitted. 

This partition function allows the calculation of the 
contribution of the fluctuations to various measurable 
quantities. For example, the contribution to the specific 
heat (per unit volume) is given in d dimensions by 



'-■a 



d'^lna{ci,i^,T) 



a/32 



Using Eq. ([24]), with t = \n{T/T^), we find 



Ofl — 



m q 



1 1 



Here, 



Cif ) - A 



q'-'dq 



1 



a(q, m,f;)2 a(q, m,f) 



(28) 



(29) 



(30) 



where A is the wave-length cutoff. 

In contrast to the 'standard' GL theory, the func- 
tion a in Eq. ([24)) grows only logarithmically at large 
|m| and |q|. Therefore, both terms in Eq. ([30)) di- 
verge unless we impose upper cutoffs on the frequen- 
cies and on the momenta. As we discuss below, some 
physical results may be affected by the resulting cutoff 
dependence. Equation (l30l) already contains the cutoff 
|q| < A. In the diffusive dirty limit, one certainly re- 
quires that |q| < hence A < l/l. As discussed in 
connection with Eq. ([T5t . we require \u\ < < l/r, 
and therefore \m\ < M = uj / {2ttT^) . Since we work 
in a regime where <^ lojj [i.e. small A, see Eq. 
([23)) ]. M is rather large. Indeed, Ref. |3 proposes us- 
ing A = l/l and M — lu j~, / {2TrT^) (but then abandons 
the quantum fluctuations altogether). However, since 
the phonon-mediated attraction arises only for momenta 
within ujjj/vp from the Fermi momentum, and since 
I = VpT and iOpi <C 1/t, one might argue that we should 



use A = Toia\{l / 1 , Lj Pi / V p} — Upi/vp. 



115 



It should be noted that in principle we has to in- 
clude also fluctuations with larger wave-lengths, up to 
the Fermi momentum pp. It has been known for a long 
time that the short-range correlations determine the crit- 
ical behavior of the internal energy. [3] In the range 
K <l/l < q <pp\.e. on distances shorter than the mean 
free path, the superconducting fluctuations are those of 
a clean superconductor. Decreasing I then yields the 



crossover from the clean to the dirty behavior. We are 
not aware of such a full calculation, and it certainly goes 
beyond the scope of the present paper. It is also not clear 
yet how to deal with the fluctuations in the intermedi- 
ate range A = ujpj/vp < q < l/l. As we show below, 
our main result is not very sensitive to the choice of A. 
In any case, our calculation gives a better description of 
the 'extremely' dirty supercunductor, when A and l/l 
approach pp . 

It is now convenient to switch to dimensionless quan- 
tities: Q = q/A, 5 = AZ and 

7 = DA^/{4:nT^) = 5"^ / {AirdT^T) 
^ ML0oT/{2d) ^ M{uja/Ep)ppl/{Ad) . (31) 

In total, the problem is described by three dimensionless 
numbers, namely 6, 7 and M, which depend on w^/T^, 
ujpi/Ep and ppl. With these, we write 



a(Q,m,t) + * 



1 



-7Q' 





1- 


- * 






.2. 



and 



a(Q,m,i)2 a{Q,m,t) 



(32) 



(33) 



Note that Cm'' depends only on t and on 7. 

With upper cutoffs, all the integrals over Q in Eq. 
(1331) remain finite, except for the 'static' terms with m = 
0. We therefore start with a detailed discussion of the 
'static' term, Cjf'^ The function a(Q,0,t) is smallest for 
small \t\ and |Q|. Defining Ql = O.OOI/7 (the prefactor 
0.001 is arbitrary, chosen so that 7Q2 ^ 1), we now 
divide the integration over Q into two regimes. In the 
first, for < IQI < Qi, we use the GL-like expansion 

a(Q,0,t) «t-HcQ2 , 



(34) 



where 



c = *'[l/2]7 = 71^7/2 , 



(35) 

^'(z) being the derivative of ^(z). Within this approx- 
imation, ^ = y/cJi/A = V T^S.ol/i&dt), is the correlation 
length associated with the fiuctuations of the Gaussian 
'static' dirty mode. Here, = Vp/T^ is the (T = 0) 
coherence length of the pure superconductor, and we as- 
sume ^0 ^ L 



We next write C, 



id) 



id) 



0,0 



id) 
0,2 



id) 
0,3 ■ 



Here and below, the subscript c stands for 



'correction'. The first term in C'-'^\ 



rid) _ 



Q'^-^dQ 

(i + cQ2)2 



(36) 



represents the leading singular contribution. Here, A — 
IJd^'"^ and 

j x'^-^dx/{l+x^f ^{2-d)Ticsc{dTT/2)/A 
Jo 

(37) 



5 



is equal to 0.5 at d = 2 and to 7r/4 for d = 1, 3. The 
correction terms include 



id) 
0,1 



Qt 



d-4 



Qi 



(t + cQ2)2 (4-d)c^ 



+ , (38) 



^0,2 — 



^a(Q, 0,^)2 yQ^a(Q,0,0) 



0{t) , (39) 



and 



C, 



0,3 



where C^;!,^, = Q'^-^dQ/it + cQ^). For d = 3, ah 
the integrals in Cq ^ converge even at i = 0. For d < 
2, Cl^% diverges at t = 0: C^^gl = ln[cQl/t]/{2c) and 



Q'^-^dQ 
Q, a(Q,0,0) 



0(t) 



(40) 



= 7r/[2(rf)i/2]-l/(cQi). This term adds a negative 

singular correction to Cq'''\ causing an increase in Cq^^ 
at small \t\. For reasons explained below, we calculate 
this term for t = = t(J ,.^^^iJ2'^/^^-'^'> [see Eqs. dB]) 
and (|46p below]. This value of t is at the border of the 
Ginzburg region, where we need to evaluate Cq'^ ■ For the 
range of parameters used below, C*^2)o,3' turns out to be 
negligible compared to the total correction term. Both 



Cq^I and Cq^^ decrease fast from large positive values as 
7 increases towards 0.5, and then decrease much more 
slowly as 7 increases above 0.5. For fixed ujjj/Ep and 
Ppl, 7 is proportional to M ~ uj^/T^ [see Eq. ([3T|) ]. so 
that it increases as decreases. 

We now turn to the quantum terms, C„(7) with m ^ 
0. All these terms have finite non-zero values at i = 0. 
Furthermore, since *[(1 + |m|)/2] - *[l/2] > 1 for aU 
TO 7^ 0, the integrand in is always negative, hence 



.(2) 



C^^ < 0. At 7 = 0, one has C, 



(d) 



{!/(*[(! +|m|)/2] 



^'[1/2]) - 1/(«'[(1 + |to|)/2] - *[l/2])2}/d, which stays 
between 0.2/d and 0.25/d for 1 < m < 10 and then 
decreases very slowly as to increases. As 7 increases, 
Crn^ decreases slowly. In any case, the sum over to grows 
as the cutoff M grows (i.e. as decreases), and adds to 
the correction term in the specific heat. 



At small T^, M is large and (unlike Cq ^ ) the total cor 



M 

25: 

m— 1 



(41) 



rection term C, 

HI 

Writing 



(d) 



increases with decreasing T^, see Fig. 



id) 



5rfA'^[At(''-4)/2 _ fJ(d)^^ 



(42) 



see Eq. (|36p. the correction Cc"^^ becomes important as 
one moves away from and t increases. Equation 



400 



200 



0.02 



0.04 




FIG. 1: The correction term ci"^^ in the Gaussian specific 
heat (thick line) and the discontinuity in the Landau specific 
heat AC'''V('S'dA'') (thin line) at d = 3 (top) and at d = 2 
(bottom), for cj^ = 400K, Ep = lO^K and ppl = 10. For 



CF' we used t = to ^ A/{TTppl) = 0.012. 



shows that the calculated C^'^ becomes negative when 
t <tQ, where 



A \2/(4-d) 



a 



a 



(d) 



T 



■'D 



D' 



d/2- 



a 



id) \ TT 



8d 2E„ 1 



d/2- 



2/{4~d) 



2/(4-d) 



(43) 



This threshold decreases rather fast as decreases. For 



d < 2, cji^^ contains the dependent term Cq"'^,, and 
one has to solve for iteratively. Since the specific 
heat should always remain positive, we conclude that the 
Gaussian approximation becomes questionable at i < ig. 
In that regime one should add higher order terms to the 
free energy. As we see below, one never reaches this 'for- 
bidden' regime, since we always have Iq > ip- 



id) 



THE GINZBURG REGION 

The fluctuations become dominant when their contri- 
bution to the specific heat Cff\t) becomes compara- 
ble to or larger than the mean-field discontinuity AC'-'*'. 
5] Comparing Eq. with AC'-'*-' gives the so-called 



6 



Ginzburg criterion, 



(44) 



Since we find that Cc''^ is mostly positive, it causes a de- 
crease in Iq, which becomes more and more significant 
as more and more quantum fiuctuations are added (i.e. 
at lower T^) (although it also contains 'classical' correc- 



tions from C| 



0,c , 



Counterintuitively, quantum fluctua- 



tions may reduce the Ginzburg regime! 

Using Eq. ((26|) . together with A = uj^j/vp, we find 



5 — TT- 



(45) 



which is linear in T^, with a large slope. 

We now consider the dependence of ig. on the tran- 
sition temperature T^. Figure 1 shows the dependence 
of both terms in the denominator of Eq. , Cc'*' and 
AC'^'^yiS^A'^), on T^. At high T^, the correction term 
is small, and one may neglect is and obtain the 'usual' 
'static' Ginzburg criterion static^ below. However, 
at low the correction term becomes large, and it can- 
not be ignored. The two curves in Fig. 1 intercept at 
crossover temperature . Although both and Ep 
are much larger than T^, usually is much smaller than 
Ep, and therefore the slope h{2Ep/ujp,Y Eq. (|i5|) be- 
comes quite large at d = 3, yielding a relatively small 
Ty^. For example, if we use Ep = lO^K, LOp, = 400K and 
Ppl = 10, the resulting value for the crossover tempera- 
ture is ~ 0.024K (see Fig. 1). However, the crossover 
temperature Ty^ becomes significantly larger at lower di- 
mensions: At d = 2 we find « IK. 

Since the cutoff on the momenta enters only via the 
dimensionless parameter 7, and since varies slowly 
with 7 at large 7, it turns out that the results for 
are not very sensitive to the value of the cutoff A. This 
is true for all ppl > 1 in d = 3, and for ppl > 2 at 
d = 2. Therefore, the questions raised above, before Eq. 
(PT|) . may not be too severe. For d = 1 the equations 
yield ^ 3K, but they also yield t^ ^ 200, which is 
certainly beyond the range of the approximations used 
above. However, the latter value is consistent with Ref. 
d, which found very large effects of the quantum fluctu- 
ations in one dimension. 

For Tj, ^ Ty we can neglect the second term in the 
square brackets in Eq. and reproduce the 'usual' 

static Ginzburg criterion. 



t. 



G, static 



2/(4-d) 



TT / 



E, 



(d~2)/2 



ippir/^ 



2/(4~d) 



(46) 

as found in many earlier papers. [T^ Interestingly, 
increases with at d = 3, does not depend on 



at d = 2 and deacreases with at d = 1. Indeed, this 
approximate expression has been adopted in most of the 
literature. 5] However, as approaches the square 
brackets in Eq. ([^^ increase, and decreases relative 
to Eq. At = we have t^ = t^^staUc/^^^^'-"^ ^ 

as used in the plot of Cc in Fig. 1. Eventually, Eq. p6)) 
is no longer valid for < Ty,. 

For < Ty, the square brackets in Eq. (|44l) are dom- 
inated by d/\ This means that the difference between 
the two terms on the right hand side of Eq. (j42|) . which 
IS of order AC^'^), is very small. As a result, the solution 
to the equation Cf^'' — AC'^'^'^ becomes very close to 
tQ, where the approach must be abandoned. Since we 
still have tQ >Ty, and we must stay above ig, we never 
encounter the regime t < t^. However, for < the 
higher order terms in Cg"*' may be important even above 
the calculated i^.. 



The 'static' t 



G^static 



and corrected (lower) t^ Ginzburg 
criteria are plotted in Fig. 2. Unfortunately, for the 
parameters used above the value of at the crossover 
point Ty for d = 3 becomes of order 10~^°, which is not 
realistic experimentally. In contrast, at d = 2 we find 
tg. « 0.012, which is quite reasonable. These numbers 
are proportional to l/{ppl)'^^^^~'^\ so they increase with 
increasing disorder. 



2. X 10"' 



i.xiO"'' 




FIG. 2: The 'static' (upper curve) and the corrected (lower 
curve) Ginzburg ranges versus the transition temperature, at 
d = 3 (top) and d = 2 (bottom), for the same parameters as 
in Fig. 1. The vertical lines show T^. 
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SUMMARY AND DISCUSSION 

Our main result concerns the crossover temperature 
Ty^ , at which the (classical and quantum) corrections 
to the leading divergent Gaussian specific heat become 
equal to the Landau mean-field discontinuity in the spe- 
cific heat. When the superconducting transition tem- 
perature is larger than , the Gaussian fluctuations 
(and especially the quantum ones) cause a (surprising) 
decrease in the Ginzburg region, by a factor which can 
be as large as 2^/^'*"''^ [see Eq. (j44|) ]. However, when 
is smaller than , the Landau mean-field discontinuity 
becomes smaller than the contribution from the fluctu- 
ations, so that the Ginzburg criterion implies an almost 
vanishing Gaussian specific heat. Therefore, one should 
no longer use the specific heat to deduce the Ginzburg 
criterion. 

In three dimensions we estimate « 0.025K, which 
is rather low. However, increases at lower dimensions, 
which casts doubts on the use of the 'standard' Ginzburg 
criterion for low transition temperature superconductors 
in those dimensions. 

It should be noted that we derived the Ginzburg crite- 
rion above T^, based on the specific heat. Other deriva- 
tions of the Ginzburg criterion are mainly below T^,, e.g. 
comparing the fluctuations in the order parameter to its 
average value, [l5l | or comparing the ordering free energy 
in a coherence volume to k^T. 17 1 In Ref. [^we compared 
the persistent current due to the Gaussian fluctuations 
to that generated by the quartic term in the GL free en- 
ergy. All of these criteria give the same scaling with the 
basic physical parameters of the problem, but yield dif- 
ferent prefactors. Future work should consider the effects 
of quantum fluctuations on these other criteria. 

It should also be noted that similar correction terms 
arise in all phase transitions (even without the quantum 
fluctuations). We hope that our paper will stimulate 
more discussion of such corrections in other cases. 

This paper is dedicated to the memory of our dear 
friend and colleague Yehoshua Levinson. We have col- 
laborated with him on a variety of physics problems, and 
we have all learned from him a lot. He is fondly re- 
membered. This work was supported by the US-Israel 
Binational Science Foundation (BSF), by the Israel Sci- 
ence Foundation (ISF) and by its Converging Technolo- 



gies Program. OEW and AA also acknowledge the sup- 
port of the Albert Einstein Minerva Center for Theo- 
retical Physics, Weizmann Institute of Science and the 
hospitality of the Korean Institute for Advanced Study. 
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